Maximally multipartite entangled states 
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We introduce the notion of maximally multipartite entangled states of n qubits as a generalization 
of the bipartite case. These pure states have a bipartite entanglement that does not depend on 
the bipartition and is maximal for all possible bipartitions. They are solutions of a minimization 
problem. Examples for small n are investigated, both analytically and numerically. 
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The characterization of multipartite entanglement is 
no simple matter. The bipartite entanglement of a com- 
posed quantum system [l| can be consistently defined 
and quantified in terms of the entropy of entanglement or 
some physically equivalent quantity. On the other hand, 
there is no unique way of characterizing the multipartite 
entanglement. Different definitions often do not agree 
with each other, because they adopt different strategies, 
focus on different aspects and capture different features of 
this quantum phenomenon 0, [3,|3|- There is a profoimd 
reason behind this problem: the number of measures (i.e. 
real numbers) needed to quantify multipartite entangle- 
ment grows exponentially with the size of the system 
(e.g., the number of qubits). Therefore, the definition 
of appropriate entanglement measures, able to summa- 
rize the most salient global features of entanglement, can 
be very difficult and their evaluation bear serious com- 
putational difficulties. This difficulty is a characteristic 
trait of complexity Q and entanglement is no exception 
in this respect 0, Q ■ we shall introduce here multipartite 
entangled states that bear the symptoms of frustration. 

The aim of this Letter is to explore "maximally" mul- 
tipartite entangled states of n qubits. These states, to be 
precisely defined later, are maximally (bipartite) entan- 
gled for all possible bipartitions. The focus is therefore 
on the global, partition-independent, features of entan- 
glement. We will consider only pure states, the extension 
to mixed states being not straightforward, due to well- 
known phenomena such as bound entanglement Q . 

We consider an ensemble S" = {1, 2, . . . , n} of n qubits 
in the Hilbert space H = (C^)®", whose state is 



E 



k^r = 1, (1) 



where k = (fci)igs, with ki £ Z2 ~ {0, 1}, and 
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(2) 



In order to analyze the multipartite features of the en- 
tanglement shared by the qubits, we proceed as follows. 
Consider a bipartition {A, A) of the system, made up of 
riA and qubits, respectively, where ^ C is a subset 
of UA elements, A = S\A its complement, ua + riA — n 
and we will stipulate riA < tia with no loss of gener- 
ality. The total Hilbert space is accordingly factorized 
into H = Ha (X) Ha, with Ha = (S'ieA dimensions 
Na = 2"-* and Na = 2"-4, respectively [NaNa = N). 
As a measure of the bipartite entanglement between the 
two subsets, we consider the purity of subsystem A 



Tr A Pa, pa = Tr^ |V'> (V-L 



(3) 



Trx being the partial trace over subsystem X. Wc notice 
that TTA — T^A ^^'^ 



1/Na < < 1. 



(4) 



State ([T|) can be written accordingly to the bipartition 

{A, A) as 



(5) 



where \k^) = {g)-g^ \ki) e Ha- By plugging Eq. (O into 
Eq. ^ we obtain 

Notice that for a given bipartition it is very easy to 
saturate the lower bound 1/Na of (O: one looks for 
those maximally bipartite entagled states that yield a to- 
tally mixed state pA = 1/Na- We will generalize the 
above property by requiring maximal possible mixedness 
for each subsystem A G S, given the constraint that the 
total system be in a pure state. A state endowed with 
this property will be called a maximally multipartite en- 
tangled state (MMES). 
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In the most favorable case this means that every sub- 
system A composed of ua < n/2 qubits is in a totally 
mixed state pA = and thus tta = I/^a (recall 

that tta = TT^, thus when ha > n/2 the above require- 
ment applies to A). In fact, it is sufficient to look at 
maximal subsystems of size ha = [n/2] {[x] =integer part 
of x), because the density matrix of every smaller part 
B C A would automatically be of the sought form, = 
TigriAPA = t/Ns- Therefore a perfect MMES _would 
be maximally entangled for every bipartition (A, A) and 
would be characterized by tta = ^/Na for all balanced 
bipartitions. Observe that the requirement of maximal 
mixedness for {A, A), tta = 1/Na, and the analogous 
requirement for a different balanced bipartition {B,B) , 
with B ^ A, might not be compatible with each other, so 
that, at variance with the bipartite case, perfect MMES 
do not necessarily exist. 

We define a MMES as a minimizer of what we shall 
call the potential of multipartite entanglement 




where nA = ['^/2]. The above quantity is related to the 
(average) linear entropy Sl = jvl^i (1 ~7'"me) introduced 
in [1] , that extends ideas put forward in ^ • See also 
0. The quantity ttme measures the average bipartite 
entanglement over all possible balanced bipartition and 
thus inherits property ([4]), i.e. 

1/Na < TTME < 1. (8) 

The upper bound ttme = 1 is attained by the fully fac- 
torizcd states, Zk = Oies '^fe ' "^^^^ lo^oP + Ic^iP = 1- 
On the other hand, the lower bound ttme = ^/Na, if 
attained, would correspond to a perfect MMES, maxi- 
mally entangled for every bipartition. However, it can 
happen that the requirements of maximal mixedness for 
different bipartitions compete with each other. In such 
a case, the system is frustrated and the minimum of the 
potential (O is strictly larger than the lower bound in 
([8|), i.e. miuTTME > ^/Na- Since in such situation it may 
happen that different bipartitions yield different values 
of TTA, our strategy will be to seek those states among 
the minimizcrs that have the smallest variance. 

This quest can be recast as an optimization problem: 
search for the minimum of the cost function 

7i'ME(A) = TTME + ActmE, (9) 

where A > is a Lagrange multiplier and 

c^ME = f ^ ] E ('^^-'^me)^ (10) 

^ ^ \A\=nA 

is the variance of tta over all balanced bipartitions. No- 
tice that the introduction of A enables one to look for 



a compromise between the minimal purity ttme (max- 
imal average entanglement) and the minimal standard 
deviation cme of the distribution (maximally distributed 
entanglement). In general, the solution of this optimiza- 
tion problem completely defines a class of states with the 
maximal possible entanglement (minimum purity), that 
is also well distributed, being as insensitive as possible to 
the particular choice of the bipartition. 

If A ^ 1 the minimization process will yield a very 
well peaked distribution of tta around its average: entan- 
glement will be uniformly distributed, but this does not 
necessarily provide a MMES; for example, a completely 
separable state has a (vanishing) entanglement that is in- 
sensitive to the change of the bipartition {tta = 1 for all 
bipartitions). More interesting is the case A = 0; indeed, 
a solution that minimizes the cost function ■?tme(0) will 
have an entanglement distribution centered on the min- 
imum of the potential ttme- Therefore, if this minimum 
saturates the lower bound in the width ctme must 
vanish. This would be our desideratum. However, it is 
known that, for ti > 8, perfect MMES do not exist 
The general problem is therefore complicated. As a first 
step we set A = in Eq. ([9|), and focus on the minimiza- 
tion of the potential ttme = 7rME(0). We shall tackle this 
problem both analytically and numerically. By plugging 
([6|) into ^ one gets after some combinatorics 

TTME = E Mki,k2;k3,k4)zk^Zk2Zk3Zki, (H) 

with 

A{ki,k2;li,l2) ^ r) Y^kf.it^i'i^ii^K^it^ki^i^ 

\A\=nA 

= g{ki ®liV k2® h, fci © ^2 V k2 © ^i), 

(;;J'w("„-J°i J"'), (12) 

where a©6 = {ai+bi mod 2)ig5 is the XOR operation, aV 
6 = {ai + bi — aibi)i^s the OR operation, aAb = {aibi)i^s 
the AND operation, and \a\ = X^ies^i- Equations (fTTj) - 
p2)) yield a closed expression for the average purity, that 
is amenable to analytic and numerical investigation. 

In order to further simplify the problem, in the follow- 
ing discussion we will replace = \zk\ in Eq. ([T]) with 
its mean value 1/\/N, and focus on the states 

Plugging Eq. ([T3)l into Eq. ([6|) we find 

iV^ + A^^ - 1 

TTA = 

N 

+ j;p E (^^m ~ ^^'rn + m' ' V'L') '(14) 

l^l' ^in-^m' 
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where (^sf^ = V'p-i(i.m), P being a permutation such that 
A = {p(l),p(2), . . . ,p(nA)}- This is an interesting for- 
mula, that is worth discussing in detail: i) first of all, 
if = X^ies'^'fc i'^ P^ - one obtains separable states, 
that yield the maximum possible value tta = 1 for all 
bipartitions: indeed the N{Na — ^ 1)/4 cosines in 

the summation in (jl4[) are all 1; ii) the first addendum in 
the right-hand side corresponds to the average entangle- 
ment of typical states 0, ■ Thus, the combination of 
phases in the summation can increase or reduce the value 
of the purity with respect to the typical one (at fixed bi- 
partition); iii) in order to get a lower value of purity, one 
should look for combinations of angles that tend to yield 
negative cosines. On the other hand, it is also clear that 
for n > 2 not all cosines can be —1, as this would yield a 
purity smaller than 1/Na', iv) at fixed bipartition (A, A) 
it is always possible to find combinations of cosines that 
saturate the lower bound tta ~ 1/Na- However, when 
plugged into (|lip , the requirement that this lower bound 
be saturated for every A might not be satisfiable. This 
problem is in general n-dependent. 

We now explicitly look at the simplest examples. For 
two qubits we only have one bipartition and the potential 
of multipartite entanglement reduces to 



.(2) _ 



1 



'ME 



= -: + -:Cos{ifa~'^i-'P2 + '^3), (15) 



where the indices of the phases are expressed again in 
terms of k, in decimal notation. The minimization of 
the potential consists in solving the equation ttJ^^e = 1/2 
(remember that Na = Na = 2). It is straightforward to 
obtain (po — fi — ^P2 + fs = tt, which yields the MMES 

In this degenerate case, multipartite entanglement coin- 
cides with bipartite entanglement, and this state is obvi- 
ously equivalent, up to local operations, to a Bell state. 
For three qubits one must look for the solutions of 



^(3) 
^ME 



1/2, where 



"me ~ g 



^ [cOs(Vp(0) + ^p{7) - fpil) ~ fp{6)) 



+ C0s((^p(2) -t- (Pp(5) - ipp(4.) - (Pp{3)) 

+ 2cos(v3p(o) + ipp(3) - <y9p(i) - (^p(2)) 

-f- 2 C0s((^p(7) + ^p(4) - (/?p(e) - </?p(5))] , (17) 

where the sum is over the 3 cyclic permutations (of the 
bits). A class of solutions is 

IV^s) = l(e*'^«|0) +e*'^i|l) +e*'^^|2) -e'(-'^»+'^i+'^^)|3) 



|4) 



|5) 



|6) 



and form a 5-dimensional submanifold. This state in- 
cludes the GHZ state [llj as a particular case and shares 
the same properties of the GHZ state for what concerns 
concurrence and one-tangle. Note that Trf^^^ contains 12 
cosines with different arguments, 6 of which are counted 
twice. The solution proposed corresponds to 2 cosines = 
1, 4 cosines = — /?, 4x2 cosines = — 1 and 2x2 cosines 
= /3, with /? = cos((^o — 'P2 — 'Pa + Ve)^ thai sum up to 
2 — 4/3 — 8 -I- 4/3 = —6. In fact, there are 3 families of so- 
lutions of the form (jl3p . corresponding to the 3 subman- 
ifolds Mp = {(^j|v3p(o) 4- (pp(7) - fp{i) - <fp{6) = 0, (pp{2) + 

fp{5) - </'p(4) - Vp{3) = 0, (^p(o) + <Pp(3) - fpil) - <fp(2) = Tt} 

with p a cyclic permutation. All three classes yield the 

(3) 

same pattern of cosines in ttj^^ (/3 being given by the 
corresponding permutation). 

For a number of qubits larger than 3, we turned to a 
numerical approach: we generated a typical state of the 
form p3p and numerically tackled the minimization prob- 
lem through different kinds of iterative algorithms (for a 
review of numerical techniques and their implementation, 
We first used deterministic algorithms. In gen- 
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oii-'Po + Vi+'Pe) 



|7» 



(18) 



eral, we found that minimization is strongly dependent 
on the initial conditions. Therefore, we sampled a large 
number of initial states, in order to test the reliability 
of the solutions obtained. Among others, the truncated 
Newton method gave us the best results in terms of both 
reliability and speed. The use of stochastic algorithms 
gave us comparable results. In both cases the existence 
of a large number of degenerate (local and global) minima 
required an accurate analysis. 

(A) 

For n ~ 4 qubits, we numerically obtained min ttj^^ ~ 
0.333 > 1/4 with ctme — 10~^. If the requirement \zk\ = 
1/Vn in Eq. ^ [and (fT3|) ] is relaxed, one can make 
Ctme vanish. This is a first example of frustration among 
the bipartitions, that prevents the existence of a perfect 
MMES. It is curious that the requirement that purity be 
minimal for all balanced bipartitions generate conflicts 
already for n = 4 qubits. This is consistent with results 
obtained by other authors l^, 14 1. 

For n = 5 and 6, the landscape of the manifold on 
which the minimization is performed becomes compli- 
cated. Nonetheless, we found perfect MMES, namely 
solutions for Tr^g = 1/4 and Tr^g = 1/8, respectively. 
Therefore, and curiously, frustration is present for n = 4 
qubits, while it is absent for n = 5 and 6. For example, 
a 5-qubits perfect MMES is defined by Eq. ^ with the 
following set of phases 

(^fc) = (0,0,0,0,0,^,7r,0,0,^,7r, 0,0, 0,0,0, 

0, 0, TT, TT, 0, TT, 0, TT, TT, 0, TT, 0, TT, TT, 0, 0) (19) 

and lives on a 7-dimensional manifold. The distribution 
of the angles x = ip^i — ip^i, +ip^iii — ^fi^ii for a given 6-qubit 
MMES is displayed in Figure [1] We observe interesting 
features, shared by all the MMES we investigated and 
for all values of n: first of all, the distribution is sym- 
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TABLE I: Perfect MMES for different n. 
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1: Distribution of the arguments x of tiie cosines in tt^^, 



for a given 6-qubit MMES. Notice the symmetries around tt 
and 7r/2. The total number of phases is 62720. Different 
phases are fuUy resolved by the binning. 



metric around x = tt; second, there is a large number of 
instances such that cosx = — 1, partially compensated by 
the contribution of cos a; = 1; third, the distribution of 
the remaining angles is symmetric around a; = 7r/2 and 
yields a vanishing contribution to ttme- 

For n = 7 we numerically found configurations with 
miuTT^Jg ~ 0.134 > 1/8 and ctme — 10~^, which im- 
proves previous bounds 0, [l3| . By minimizing the cost 
function ^ with a non vanishing A, one can reduce ctme 
to ~ 10""^, at the expense of a higher ttme — 0.136. It 
is not clear at present if the impossibility to reach the 
absolute minimum for n = 7 is to be ascribed to the 
numerical procedure. For 7i = 8 and 9, where perfect 
MMES do not exist [3| , the convergence of the numerical 
simulations become very slow. This is a typical signature 
of frustration. The numerical results show that, in order 
to obtain a vanishing width of the distribution for frus- 
trated systems, it is necessary to increase the value of the 
average of the purity. These conclusions arc summarized 
in Table U 

A comment is now in order. Although in this article 
perfect MMES are exhibited as minimizers of the po- 
tential of entanglement, they are, by their very nature, 
independent of the method used to find them. In fact, by 
virtue of their maximal mixedness, they saturate all mea- 
sures of multipartite entanglement. This is the case, for 
example, of the global entanglement measure of Meyer 
and Wallach Q, as well as its generalizations. Indeed, 
any entanglement monotones, being functions of the par- 
tial density matrices of subsystems of qubits, attain their 
maximal values on perfect MMES. The minimization pro- 
cedure we propose is just a convenient way to construct 
them. For non-perfect MMES, due to a possible finite 
value of (Tme, different measures of entanglement can re- 
solve a part of the degeneracy of the manifold of mini- 
mizers. 

It is interesting to briefly discuss one straightfor- 
ward potential application. Consider the 5-qubit perfect 
MMES I^Z-s) defined by Eqs. dTS]) and (HH). One can easily 
prove that 



n 


perfect MMES 


2,3,5,6 
4 
7 
> 8 


exist 
do not exist 
? 

do not exist 



where a\ and erf are Pauli matrices. Therefore, the 
single- and two-qubit statistics are always "flat", but 
the measurements of the observables in (pO)) are always 
strictly correlated. Any two parties, that can be far 
apart, can therefore share a cryptographic key only if 
the other three parties agree on measuring their respec- 
tive observables in (j20p and making their results public. 
Notice that the key is shared by the two parties but is 
unknown to the other three. We call this phenomenon 
majority-agreed key distribution. 

In conclusion, we introduced a class of multipartite 
entangled states that maximize the amount and distri- 
bution of entanglement. The features of these states 
strongly depend on the number of qubits involved. In 
our numerical search, wc noticed that already for a rela- 
tively small number of qubits (n > 7), the landscape of 
the parameter space where the optimization procedure 
is performed has a complex structure with a large num- 
ber of local minima and a very slow convergence. The 
presence of frustration, due to the competition among 
different partitions (observed already for n = 4), appears 
to be a general feature of many-body systems. It prevents 
the possibility to find perfect MMES but introduces in- 
teresting perspectives. In this sense the minimization 
task is a problem that requires a careful analysis and the 
use of numerical and analytical strategies from different 
research fields. In this paper we just started to explore 
these connections. Thus, the study of maximally multi- 
partite entangled states paves the way towards a deeper 
comprehension of the complex structure of quantum cor- 
relations arising in many-body systems. 

This work is partly supported by the European Com- 
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(20) 



W. K. Wootters, Quantum Inf. and Comp., 1, 27 

(2001) ; M. B. Plenio and S. Virmani, Quant. Inf. 
Comp. 7, 1 (2007); L. Amico, R. Fazio, A. Osterloh 
and V. Vedral "Entanglement in Many-Body Systems," 
arXiv:quant-ph/0703044 (Rev. Mod. Phys., in print). 
V. Coffman, J. Kundu and W. K. Wootters, Phys. Rev. 
A 61, 052306 (2000); A. Wong and N. Christensen, Phys. 
Rev. A 63, 044301 (2001); D. Bruss, J. Math. Phys. 43, 
4237 (2002); G. Rigolin, T. R. de Oliveira and M. C. de 
Oliveira, Phys. Rev. A 74, 022314 (2006). 

D.A. Meyer and N.R. Wallach, J. Math. Phys. 43, 4273 

(2002) . 

A. J. Scott, Phys. Rev. A 69, 052330 (2004). 



5 



[5] K. R. Parthasarathy, Proc. Indian Acad. Sciences, 114, 
365 (2004). 

[6] M. Mezard, G. Parisi and M. A. Virasoro, 5pm Glass 
Theory and Beyond (World Scientific, Singapore, 1987). 

[7] V. I. Man'ko, G. Marmo, E. C. G. Sudarshan and F. 
Zaccaria, J. Phys. A: Math. Gen. 35, 7137 (2002). 

[8] P. Facchi, G. Florio and S. Pascazio, Phys. Rev. A 74, 
042331 (2006); Int. J. Quantum Inf. 5, 97 (2007). 

[9] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 
Rev. Lett. 80, 5239 (1998); C. H. Bennett, D. P. Di- 
Vincenzo, T. Mor, P. W. Shor, J. A. Smolin, and B. M. 
Terhal, Phys. Rev. Lett. 82, 5385 (1999). 
[10] E. Lubkin, J. Math. Phys. 19, 1028 (1978); S. Lloyd and 
H. Pagels, Ann. Phys., NY, 188, 186 (1988); D. N. Page, 
Phys. Rev. Lett. 71, 1291 (1993); K. Zyczkowski and H.- 



J. Sommers, J. Phys. A 34, 7111 (2001); A. J. Scott and 

C. M. Caves, J. Phys. A; Math. Gen. 36, 9553 (2003). 
[11] D.M. Greenberger, M. Horne and A. Zeilinger, Am. J. 

Phys. 58, 1131 (1990). 
[12] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. 

Flannery, Numerical Recipes in Fortran 77, Numerical 

Recipes in C (Cambridge University Press, 1992). 
[13] A. Higuchi and A. Sudbery, Phys. Lett. A273, 213 

(2000); S. Brierley and A. Higuchi, J. Phys. A: Math. 

Theor. 40, 8455 (2007). 
[14] A. Borras, A.R. Plastino, J. Batle, C. Zander, M. Casas, 

and A. Plastino, J. Phys. A: Math. Theor. 40, 13407 

(2007). 



